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ABSTRACT 



The objectives of the project described in this 
thesis were to ; 1) provide the documentation that is 
needed for a Naval Postgraduate School student to use 
the general purpose finite element computer program 
called SAP IV, and 2) to make available, and prepare 
the users manual for, a pre-and postprocessor program 
called SOBROUTINE PSAP. This subroutine, which was 
developed at the NASA Langley Research Center, has 
been modified to specifically plot the finite element 
model gecmetry for SAP IV models and to postprocess 
displacement data for those models on the NPS Calcomp 
Model 765 Plotter, The input and output for SAP IV and 
SUBRODTINE PSAP are discussed in detail. The codes 
have been used successfully in AE 4102, Flight Vehicle 
Structural Analysis II. 



4 



TABLE OF CONTENTS 



I. INTRODUCTION 8 

A. STRUCTURAL ANALYSIS SOFTWARE PACKAGES 8 

3. PRE-AND POSTPROCESSORS FOR STRUCTURAL ANALYSIS 

PROGRAMS 9 

C. THESIS MOTIVATION 9 

II. GUIDE TO THE USE OF SAP IV AT NPS 11 

A. DESCRIPTION OF SAP IV 11 

B. COMPUTER CARD DECK PREPARATION AT NPS FOR 

SAP IV 12 

C. HELPFUL POINTS ON DATA DECK PREPARATION 15 

D. ALTERATION OF SAP IV AT NPS 16 

III. GUIDE TO THE USE OF SUBROUTINE PSAP, A PRE-AND 

POSTPROCESSOR FOR SAP IV 20 

A. SUBROUTINE DESCRIPTION 20 

B. SUBROUTINE PSAP OSS 21 

1. General Set-up of Input Deck 21 

2. Job Control Cards 22 

3. Main Program 25 

4. Title Card 25 

5. Namelist OPTION 26 

6. Geometry Input Data Deck 28 

7. Case Indentif ication Card 31 

8. Deck of Displacement Data to be Plotted.. 33 

9. Namelist PICT 34 

C. METHOD FOR ALTERING SUBROUTINE PSAP 38 

D. SIGNIFICANT ASPECTS OF SUBROUTINE PSAP 39 

1. Exploded Plots 39 

2. Portions of Models 39 

3. Specification of View 40 

4. Possible Sources of Errors 40 



5 



IV. CONCLUSIONS AND RECOMMENDATIONS 42 

Appendix A: EXAMPLE PROBLEM USING SAP IV 44 

Appendix B; EXAMPLE PROBLEM USING SUBROUTINE PSA? 46 

Appendix C: SAMPLE SUBROUTINE PSAP OUTPUTS FOR A FINITE 

ELEMENT MODEL OF A WING 53 

Appendix D: SUBROUTINE PSAP LISTING 59 

LIST OF FIGURES 60 

LIST OF REFERENCES 89 

INITIAL DISTRIBUTION LIST 90 



6 



ACKNOWLEDGEMENT 



I wish to acknowledge the help that I have received in 
the course of my research. 

I am particularly indebted to Professor R. E. Ball, who 
provided me with the guidance and aid necessary to complete 
the project. 

I am also grateful to Professor Gilles Cantin, who 
provided me with valuable information in starting and 
carrying through my work. 

I also wish to give credit to the NPS programming 
consultants, and in particular Sharon D. Raney, for the fine 
assistance I received in debugging my computer runs. 

My gratitude is also heartfully extended to Richard 
Citerley for providing me with the computer card decks 
necessary to begin my project. 

I especially want to express my sincere appreciation to 
my wife, Karen, for her patience and perseverance with the 
typing of this thesis. 



7 



I. 



INT RO DUCT ION 



A. STRUCTURAL ANALYSIS SOFTWARE PACKAGES 



The past decade has seen great strides in the field of 
computer software packages that have been developed for use 
in structural analysis. A great deal of time and effort has 
been expended in the development and distribution of these 
packages. Today there are structural analysis computer 
programs for almost every conceivable structure an engineer 
could desire to analyze. The types of analyses performed by 
a given package vary widely and may include features for 
linear or nonlinear materials, static analysis, dynamic 
analysis, buckling analysis, and nonlinear dynamic analysis, 
to name only a few. These programs can be grouped into two 
major categories - special purpose or general purpose. 

In order to expose the students at the Naval 
Postgraduate School to the use of a general purpose 
structural analysis program, as well as to provide the 
capability of using such a program in research work, the 
finite element structural analysis program SAP IV (Reference 
1) was acquired and made operational at the Naval 
Postgraduate School by Professor Gilles Cantin of the 
Mechanical Engineering Department. SAP IV can perform 
linear static and dynamic analyses on one-, two-, and 
three-dimensional structures. 
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B. 



PRE-AND POSTPROCESSORS FOR STRUCTURAL ANALYSIS PROGRAMS 



After the development of a large number of structural 
analysis programs, users began to recognize that a disparity 
existed between efficient general purpose structural 
analysis programs and optimum utilization of these programs. 
Many of today's software structural analysis packages, such 
as SAP IV, reguire the user to prepare and reduce tremendous 
amounts of data. The need existed for some aids in 
processing and reducing these large quantities of data. 
Consequently, there have been many pre- and postprocessors 
developed for a specific use as well as for a general use 
basis during the past several years. The value of a given 
processing package lies in its ability to aid the user in 
preparing his model geometry, in his data checks, and in 
processing the output in an easily understood fashion. One 
of the most effective means utilized in preparing or 
reducing data is through the use of visual displays, whether 
they be designed primarily for graphic presentations, such 
as the Strcmberg-Carlson, or for paper plots, such as 
Calcomp. This ability to visually display input and output 
data is a highly valuable tool for the structural analyst. 



C. THESIS MOTIVATION 



The desire to have a pre- and postprocessor that could 
be used in conjunction with the general purpose structural 
analysis program SAP IV prompted the acquisition and 
implementation of a general use plotting package by this 
author. After researching the possible options, an existing 
program was obtained from Anamet Laboratories, San Carlos, 
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California. The program was originally developed at the 
Langley Research Center, Hampton, Virginia by Gary L. Giles 
for use with modern digital computers. The program, details 
of which can be found in Reference 2, generates oblique 
orthographic projections of three-dimensional finite element 
models and is distinguished by its provisions for 
generality, ease of use, different display options, and 
computational speed. The computer code was written for use 
on CDC 6000 series machines and had to be modified somewhat 
for use on the NFS 360/67. The modified version of the 
program is now available for use with the NFS Calcomp 
Flotter Model 765. 

The primary purpose of this thesis is to provide the 
necessary documentation in order that students enrolled at 
the Naval Fcstgraduate School, and specifically in the 
courses AE 4101, 4102 (Flight Vehicle Structural Analysis I 
and II respectively ) , can, with a minimum of difficulty, 
effectively utilize SAF IV and its now- functional pre-and 
postprocessor FSAF. The remainder of this thesis is broken 
down into two major subdivisions, 

1) II. GUIDE TO THE USE OF SAP IV AT THE N AVAL 
POSTG RADUATE SCHOOL 

2) III. GUI^ TO THE U^ OF SUBROUNT INB PSAP, A 
PR E-AND POSTPROCES SOR FOR SAP IV 

Appendices A and B of this work provide detailed examples 
for the input preparation and output reduction of data using 
both SAP IV and SUBROUTINE PSAP. 
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II. GOIDE TO THE USJ OF SAP IV ^ NPS 



A. DESCEIPIldN OF SAP IV 



SAP IV is a general parpose structural analysis digital 
computer program that can provide a finite element solution 
for both the static and dynamic analysis of linear 
structural systems. A. detailed user's manual is contained 
in Reference 1. The program has the capacity to analyze very 
large three-dimensional systems, as well as small systems, 
with no loss in efficiency. SAP IV, which is coded in 
FORTRAN IV , is a very flexible program and can be 
considered a very efficient aid to the analyst. The purpose 
of this section is to provide the necessary additional 
documentation, above that provided in Reference 1, for a 
student at the Naval Postgraduate School to make use of the 
program. 

The methods of analysis and the construction of the 
program are not included in this section, but can be found 
in Reference 1. The program contains nine finite elements of 
the following types: 

(a) three-dimensional truss element, 

(b) three-dimensional beam element, 

(c) plane stress and plane strain element, 

(d) two-dimensional solid element, 

(e) three- dimensional solid element, 

(f) variable-number-nodes thick shell and 
three-dimensional element, 

(g) thin plate or thin shell element. 



(h) boundary element, 

(i) pipe element (tangent and bend) . 

There are numerous options and combinations of static and 
dynamic analysis that are available to the user of the 
program. Reference 1 provides specific details of the many 
available user options. 

B. COMPUTER CARD DECK PREPARATION AT NPS FOP. SAP IV 

Figure 1 outlines the overall computer card deck 
necessary to access and utilize SAP IV as it is currently 
operational at NPS. A complete detailed breakdown of the 
necessary IBM job control cards follows on the next page. 
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// ( STANDARD GREEN JD3 CARD ) 

//GO EXEC PGM=SAP,REGI0N=260K 

//STEPLIB OD UN IT=2321 ,VOL=SER=CEL002 ,DISP=SHR, 
// DSN=F0559.SAPLM 

//GO.FTCIFOOI DO UN IT=SY SD A , SPAC£= { CYL » ( 6, 6 ) ) , 
// 0CB=(RECFM=VS3,LRECL=1284,BLKSIZE=2572J 
//G0.FT02F001 DD UN IT = S Y SO A, SPACE= ( CY L , (6 , & ) ) , 
// DCB=(RECFM=VSB,LRECL=1284,BLKS IZE=2572J 
//GD.FT03F001, DD UN IT = S Y SO A , SPACE = ( C Y L , ( 6 ,6 ) ) , 
// DCe=(RECFM=VSB, LREC L=1284, BLKS IZE = 2572) 

//G0.FTG4F00L DD UN I T = S Y SD A , SPAC E= ( CYL » ( 6 , 6 ) ) , 
// DCB=(RECFM=VSB,LRECL=1284, BLKS IZE=2572 ) 
//GQ.FTC5F001 DO DDNAN1E = SYS IN 

//GC.FTC6F001 DD SYSOUT=A , S P ACE = ( CYL , ( 3 , 1 ) ) , 

// 0CB=(RECFM=FBA,LRECL=133,BLKSIZE=1330) 
//GQ.FTC7F001 DD UN IT=SY SDA , SPACE= ( CYL , (6 ,6 ) ) , 
// 0CB=(RECFM=VSBtLRECL=l284,BLKSIZE=2572) 
//G0.FTC8F00L DD UN I T = S Y SO A , SPAC E= ( C YL , (6 ,6 ) ) , 
// DC6=(PECFM=VSBt LRECL=1284,3LKSIZE=2572) 
//G0.FTC9F001 DD UN I T = SY SO A , SPACE = ( C YL , ( 6 , 6 ) ) , 
// DC6=(RECFM=VSBtLRECL=1284,BLKSIZE=2572) 
//GQ.FTIOFOOL DD UN IT=SY SD A , SPACE= ( C YL , ( 6 , 6 ) ) , 
// DC E= (RECFM=VSBt LREC L= 1284, BLKS I ZE=2572 ) 

//GO.FTllFOOl OD SYSaUT=8,SPACE=(TRK, (20,2)) , 

// DCB=(RECFM=FB,LRECL=80,BLKSI ZE=7200) 
//GG.SYSIN DD 

( SAP IV DATA ) 



/ ( STANDARD NPS EOF CARD ) 
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After the //GO.SYSIN DD * card, the deck of SAP IV data that 
is prepared according to the instructions in the appendices 
of Reference 1 follows. Following the SAP IV data is the 
standard NPS end-of-file card (/♦) . 



C. HELPFUL POINTS ON DATA DECK PREPARATION 



There are several possible areas where a user can err in 
his preparation of data for SAP IV. The following items are 
representative of a few common sources of error. 

1) Particular attention should be directed toward 
using the correct formats and correct card columns in key 
punching data on cards. (i.e., integer formats, right 
justified) 

2) The program has several internal data generation 
features inherent in it. Simply because data were generated 
during the program execution does not necessarily mean they 
were generated correctly. Errors in the user- prepared input 
cards used in data generation can cause severe discrepancies 
to occur during program execution. The generated data 
should be carefully checked for accuracy. 

3) In order to terminate a given problem a number of 
blank cards are necessary at the end of the input data deck. 
Specific points of interest pertaining to the problem of 
termination can be found in Reference 1 under note (1) on 
page II. 1 and in note (1)* on the top of page V.2. 

SAP IV has a wide range of options where a clever user 
can fully exploit the full capabilities of the program. The 
detailed decription of these features is found in Reference 
1; however, the following list summarizes a few of the more 
useful options and features. 
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1) Data check only mode of execution. 

2) Nodal point and element data generation 
capabilities. 

3) Five different types of analysis. 

a. static analysis 

b. eigenvalue /vector solution 

c. forced dynamic response by mode superposition 

d. response spectrum analysis 

e. direct step by step integration 

4) Automatic punched computer card output of 

displacement results. This feature was added to the program 
by this author to make possible graphic postprocessing. 

The foregoing discussion is not intended to be a 
complete diagnostic summary of SAP IV, but rather an aid to 
the student who desires to get started using the program as 
it currently exists on the NPS IBM 360/67. The above 
discussion, along with a copy of Reference 1, should provide 
a jumping off place for a novice structural analyst. A 
complete example of the job control cards and input data 
deck for a static truss analysis can be found in Appendix A. 



D. ALTERATION OF SAP IV AT NPS 



Because SAP IV is a very flexible program, it is 
possible for a user to make modifications to the basic 
program. The computer code is complex, but it is not overly 
difficult to modify parts of the program in order to satisfy 
a specific user need. A method that has proven successful 
for this author, in the modification of a SAP IV subroutine 
to provide punched output of displacement data, and to 
create a personal version of SAP IV, is outlined in this 
section. 
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The first step is to define the subroutine name that the 
user wishes to modify. A listing of the source program is 
available through Professor R. E. Ball, Department of 
Aeronautics, or Professor G. Cantin, Department of 

Mechanical Engineering. Having defined those portions of the 
program for which modification is desired, the next step is 
to obtain a punched copy of those desired routines by using 
the following format. 

// (STANDARD NPS JOB CARD) 

//SYSPRINT DD S YSOUT=A, S P ACE= (TRK , ( 1 0 , 1 ) ) 

//SYSUT1 DD DISP=SHR,aNIT=2321,VOL=SEE=CEL001, 

// DSN=F0099. SAPSR 
//SYSDT2 DD SYSOOT=B 
//SYSIN DD 

PONCH TYP0RG=P0,MAXNAflE=2 
MEMBER NAM E= (SUBROUTINE NAME) 

/* 

NOTE: SUBROUTINE NAME is the name of the SAP IV subroutine 
required . 

With the desired routine decks now in hand, the 
necessary changes can be incorporated into the subroutine 
deck and an object deck is then obtained as follows: 

// (STANDARD NPS JOB CARD) 

// EXEC FORTCD 
// FORT. SYSIN DD * 

(MODIFIED FORTRAN SOURCE DECK) 

/* 

The next and final step is to take the object deck 
obtained in the previous step and insert it into the proper 
position in the control cards that are illustrated on the 
following page. 
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)ARC NPS EOF CARD 



There are many options available in the creation of load 
module libraries and they are discussed in detail in 
Reference 4, sections II. and III. The portion of the 
preceding example control card deck that would necessitate 
modification is the //LINK.SYSLMOD card, where 
F0559 . SAPLM (SAP) should be changed according the following 
format : 

General format: Lnnnn. anyname 
where L is : 

S-student data set 
F-faculty data set 
nnnn is : 

user number assigned 
anyname is : 

any unique name assigned by the user (1-6 
characters in length with the first character alphabetic ) , 
and CEL002, should reflect an appropriate data cell with 
available space. 

The procedure outlined above will create a modified 
version of SAP IV on a chosen data cell available to the 
user for a period of 90 days. The program can now be 
executed with appropriate modifications to the //STEPLIB 
card as discussed previously in section II. B. Any questions 
concerning the creation of load module libraries or their 
execution can be answered by any of the programing 
consultants on the first floor of Ingersoll Hall. 
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III. 



guide to the use of subr outine psap, a pre-and 

POSTPROCESSOR FOR S^ 17 



A. SUBROUTINE DESCRIPTION 



SUBROUTINE PSAP is a modified version of the oblique 
orthographic projection program that is found in Appendix B 
of Reference 2. The program, originally developed for use 
at the NASA Langley Research Center, Hampton, Virginia, 
required some changes so that it could be used in 
conjunction with SAP IV and the NPS Model 765 Calcomp 
Plotter. The original version of the plotting package 
allowed for various geometry and displacement data input 
options. However, the subroutine, as currently filed, is 
constructed strictly for use with the input and output of 
SAP IV. With some slight modifications to the subroutine, it 
could be adapted to any number of different types of input 
geometry or displacement data decks. A method for 
accomplishing this is outlined in section III. C. 

The current capability for generating oblique 

orthographic projections of SAP IV finite element models is 
limited to the following types of elements: 

1) Type 1, Three-Dimensional Truss Elements 

2) Type 2, Three-Dimensional Beam Elements 

3) Type 3, Plane Stress Membrane Elements 

4) Type 4, Two-Dimensional Finite Elements 

5) Type 6, Plate and Shell Elements (Quadrilateral) . 
The undeformed topology of the finite element model, useful 
in checking input data, as well as the displaced topology 



20 



projection of the same model, can be obtained from PSAP. The 
subroutine contains many different options and permutations 
of those options, some of which are listed below. 

1) plots of models annotated with grid point numbers 
of element numbers 

2) plots of portions of models 

3) exploded plots of model sections (i.e., line 
elements coincident with the edges of triangular or 
quadrilateral elements may be difficult to single 
out; program provides a capability to separate 
elements so that their absence or presence is easily 
detectable) 

4) displacements superimposed on grid point 

coordinates of the undeformed structure 

5) displacements represented as vectors extending the 
undisplaced grid points. 

There are many more combinations of options that are 
available and will be more specifically outlined in the 
remainder of this section. Appendix B details a complete 
input-output example of a pin- jointed truss. 



B. SUBROUTINE PSAP USE 



1 . Genera l S^t^u£ of Input Deck 

In general, the correct sequence of computer cards 
required to utilize SUBROUTINE PSAP is shown schematically 
in Figure 2 and consists of eight separate major groups as 
follows; 

1) a group of JCL cards and main program to allocate 

storage 

2) a single card containing title information 

3) Namelist OPTION containing values to determine if 
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proper storage allocation is available and specifying 
various program options 

4) a geometry deck (SAP IV data deck) containing grid 
points and connectivity of the model 

5) an optional single title card used to identify the 
deck of displacement data to be plotted 

6) a single card containing the value of the total 
number of SAP IV load cases and an optional scale 
factor 

7) the deck of displacement data to be plotted 
(output of the execution of SAP IV-static analysis) 

8) Namelist PICT containing values to specify the 
type of plot desired and what information is to be 
included on the plots. 

By repeating parts of the basic input data to the program, 
different plots of the same data can be generated. 



2. Job Contro l Card s 



The sequence of cards on the following page depicts 
the necessary JCL to execute SOBROOTINE PSAP. 
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// (STANCARD GREEN JOB CARD) 



// EXEC FORTCLGP, REGION. G0=180K 
//FCRT.SYSiN CD 

( MAIN PROGRAM ) 

/ ( STANCARC NPS EOF CARD ) 

//LINK.USDD DC UNIT=2321 , V0L=SER=CEL002,0I SP=SHR , 
// 0SN=F0559. PSAPLM 
//LINK'.SYSIN CD 

INCLUDE USDD(PSAP) 

ENTRY MAIM 

/ ( STANCARC NPS EOF CARD ) 

//GO.FTlOFOOl DO UNIT=SYSDA, 

// SPAC£ = (CYLt(3,l) ) , 

// CCE=(RECFM=VS»BLKSIZE=3520) 

//GG.SYSIN DC ^ 

( DATA STARTING WITH THE TITLE CARD SE 

/ ( STANCARC NPS EOF CARO ) 



FIGURE 2) 
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^ • Main P rogram 



The main program consists of three cards. It 
allocates the proper storage and calls SUBROUTINE PSAP. An 
example of a main program for a finite element model that 
contains 400 nodes is : 

DIMENSION ZZZ (2800) , DISPD (5,3,400) 

CALL PSAP (ZZZ, 2800, DISPD, 400) 

END 

The definition of the arguments used in calling SUBROUTINE 
PSAP (ZZZ, NZ, DISPD, NON) are : 

ZZZ-blankcommon array used to store nodal coordinates 
and displacements 

NZ-length of array ZZZ (NZ is determined by multiplying 
the number of nodes in the model by seven, i.e., NZ=no. of 
nodes * 7) 

DISPD-a three dimensional working array used in 
subroutine DATA9 

NON-number of nodes in the model 
It is crucial to dimension ZZZ (7*NON) and DISPD (5, 3, NON) 
correctly in the main program. Improper dimensioning of 
these arrays can cause output errors that are not readily 
traceable. 



4 . Tit 1 e ^rd 



This card contains any desired alphanumeric 
information in columns 1 to 80. The title will appear at the 
beginning of the plots. 
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5 . Mie li^ 



Namelists are a convenient means of inputing the 
names of several parameters along with their corresponding 
values. ?or the NPS IBM 360/67, the format for using 
namelists is as follows: 

1) card 1-&name-beginning in card column 2 where name 
is the name of the subject namelist 

2) succeeding cards- beginning in card column 2, the 
names of the variables and their values separated by 

, commas 

3) final card -SEND-start ing in card column 2. 

An example of a namelist format is shown below: 

col 

2 

60PTI0N 

NNDEST = 400,NaDISP =1, 

PSIZE =8.0 
SEND 

Any or all parts of a defined namelist may be included, and 
each parameter may be specified in any order between the 
Sname card and the SEND card. 

The description of the variable names in Namelist 
OPTION and their default values are contained in Reference 
2. They .are given here to assist the user in data 
preparation. 

FORMAN ^me - D efault value - D escr ip tio n 
NNDEST - 1 The number of nodes (NON) as defined in the 

program 

NUDISP - 1 0 x-direction displacements not to be input 

1 x-direction displacements to be input 
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NVDISP - 1 
NWDISP - 1 



KGEOM - 9 



KATA - 9 



NVALOS - 0 
lEESEQ - 1 



KPLOT - 1 
XSPACE - 5.0 



0 y-direction displacements not to be input 

1 y-direction displacements to be input 

0 z-direction displacements not to be input 

1 2 -direction displacements to be input 

( NOTE: when SAP IV displacement data is to 
be used, NUDISP=N7DISP=NWDISP=1 ; for no 
displacement data NUDISP=N VDISP=NHDISP=0) 
Specifies the subroutine and corresponding 
method of input for model geometry 

1 subroutine GEOm, a user-supplied 
subroutine 

2 subroutine GE0a2, a user-supplied 
subroutine 

9 subroutine GEOM9, reads in grid points and 
element data specifically from a SAP IV data 
deck 

Specifies the subroutine and corresponding 
method of input for displacement data 

1 subroutine DATA1, a user-supplied 
subroutine 

2 subroutine DATA2, a user-supplied 
sub routine 

9 subroutine DATA9, reads a punched output 
displacement deck from execution of SAP IV 
NOT INC0RP02ATED-ALL0W DEFAULT 
Grid point numbers are stored in the program 
from 1 to the total number of grid points 

0 no internal resequencing of grid points 
necessary; they are already in ascending 
order starting with 1 

1 resequence grid points from lowest grid 
point number to highest grid point number 
Specifies the type of output device to be 
used ( ALLOW DEFAULT ) 

Space between plots in the y-direction , in 
inches ( see Figure 3 for an explanation of 
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PSIZE 



10.0 



axis orientation ) 

Paper size in x-direction, in inches ( used in 
scaling of plots to insure this dimension is 
not exceeded) 

IDCASE - 0 0 no identification card proceeds the deck 

of displacement values 

1 identification card proceeds the deck 
of displacement values 



6. Ge ometr y Input Data Deck 



This portion of the input deck contains the grid 
point locations and the element connectivity. The deck has 
one of the following forms, depending on the value of KGEOM 
in the Namelist OPTION. 

KGEO M = 9 

Calls subroutine GE0M9, which is constructed to read 
SAP IV geometry data. 

(a) When KGEOM is specified as 9, the complete input 
deck of SAP IV data, prepared according to Appendices I 
through IV of Reference 1, is placed after the &END card of 
Namelist OPTION. The portions of SAP IV data deck that 
involve load, mass, or dynamic analysis data are not part of 
the input geometry data to SUBROUTINE PSAP. Only the grid 
point locations and the element connectivity data are used 
to generate the orthographic projection of the model. 
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Figure 3 - NFS CALCOMP PLOTTER ORIENTATION 
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(b) It is possible to input only a portion of the 
finite element model for a data check. To do so, it is 
necessary to modify SAP IV element control cards (described 
in detail in Reference 1, Appendix IV) for the desired 
element types to reflect the portion of the element 
connectivity cards that will be input. For example, if only 
connectivity for element numbers 15 through 50 of element 
type 1 (truss elements) of a SAP IV data deck are available 
for input, it is necessary to alter the element control card 
for the truss elements. The field on the SAP IV control card 
that defines the total number of truss elements (card 
columns 6-10) would reflect the upper bound, in this 
case-50, and card columns 65-70 would reflect the lower 
bound, in this case- 15. All nodal coordinates for the 
entire model may be input, or only those that specifically 
define the portion of rhe finite element model to be 
plotted. In either case, the nodal coordinates that relate 
to element numbers 15-50 must be specified. Unknown results 
will occur when trying to plot an element whose node points 
are not specified. The above feature is valuable in a case 
where several different people are preparing different parts 
of a large data base for a SAP IV problem and desire to 
individually check their inputs graphically for accuracy. 
KGE OM = 2 

Calls subroutine GE0M1 , which is prepared by the user 
to read geometry data from a program other than SAP IV. 

KGEOa f 2 

Calls subroutine GE0M2, which is prepared by the user 
to read geometry data from a program other than SAP IV. 

Use of KGE0M=1 or 2 requires modification of 
SUBROUTINE PSAP to fit the specific format of the user’s 
input geometry data. A method for doing this will be 
discussed in paragarph C. of this section. 
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7. C ase Inden tif ica tion Car d 



If IDCASE =0 is specified in Namelist OPTION, this 
card is omitted. The card, if present, contains any desired 
alphanumeric information in card columns 1-80 which will 
identify all displacements for a given case. For IDCASE = 1 
and SAP IV punched displacement data, a case identification 
card must appear before each Namelist PICT for every 
different load case that is plotted in addition to load case 
one. This is illustrated in Figure 4. A maximum number of 
five different load cases can be obtained from SAP IV. The 
case ID card information, if present, will appear before 
each load case's DISPLACEMENT DATA TO BE PLOTTED section in 
the printed computer output from SUBROUTINE PSAP. This 
information does not appear on any Calcomp plots. 
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Figure 4 - SAMPLE SUBROUTINE PSAP EXECUTION DECK WITH 
IDCASE=1 
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• Deck of Dis£^cement D ata to be Plotted 



This deck contains the displacement values of the 
nodal points. The deck has one of the following forms, 
depending on the value of KDATA specified in Namelist 
OPTION. 

KE ATA = 9 

Calls subroutine DAIA9, which reads SAP IV displacement 

data. 



(a) A single card, format (I5,F10.0), containing 
the number of total load cases (NTLC) that are in the SAP IV 
output displacement data deck, and a scale factor (SCAL3F) 
that is used in scaling the displacement data, must be input 
before the displacement deck. This card controls the input 
of SAP IV displacement data through subroutine DATA9. For 
example, a displacement deck from an execution of SAP IV 
that contains the maximum of five load cases would have NTLC 
= 5. The scale factor could be any value desired by the 
user, with default equal to 1. The DATA9 control card for 
the above case would be a 5 in card column five and the 
desired scale factor (SCALEF) in card columns 6-15. 

(b) The deck of displacement data obtained from SAP 
IV follows the DATA 9 control card. Since the maximum 
number of load cases that are punched out by SAP IV is five, 
the maximum NTLC is five and the actual number must be 
specified. Note that if NUDISP, NVDISP, and NWDISP are all 
specified as zeroes through default or in Namelist OPTION, 
the displacement data deck and the DATA 9 control card are 
not required. This feature enables preprocessing only of a 
given finite element model by PSAP. 
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(c) The parameter DMAGS in Namelist PICT also 
provides for magnification of displacements. (See section 
III.B.9) 

= Jl 

Calls subroutine DATA1, which is prepared by the user 
to read displacement data from a program other than SAP IV 
to be plotted. 

OATA = 5 

Calls subroutine DATA 5, which is prepared by the user 
to read displacement data from a program other than SAP IV 
to be plotted. 

Use of KDATA =1 or 5 requires that SUBROUTINE PSAP 
be modified to fit the format that the user's input 
displacement data specifically requires. A method for 
accomplishing this change will be discussed in parapgraph C. 
of this section. 



9 . list PICT 



The format for Namelist PICT is the same as that 
required for Namelist OPTION. It requires a single card, 
SPICT, followed by the specified parameter cards and the 
SEND cards, all cards beginning in card column 2. This 
namelist contains the values needed to specify the type of 
plot that is desired and what information is to be included 
on the plots. A detailed summary of Namelist PICT is 
contained in Reference 2 and is given here for user 
convenience. 

FORT R AN naSJ “ Defa u lt y_alue - Descript ion 
KHORZ - 1 Integer designating the horizontal axes of 

the viewing plane where , 2=1^, and 3 =Zq. 

(see Figure 5.) 



34 




35 



ORTHOGRAPHIC PflOJ. SHOWN IN X - Z VIEWING PLANE 



KVEET - 2 



PHI - 0.0 



THETA - 0.0 



PSI - 0.0 



NEWER - 1 



ISCALE - 1 

PLOTSZ - 10.0 
XORGN - 0.0 
YOEGN - 0.0 
PSCALE - 1.0 

NOTAT - 0 

XLHT -0.15 
KDISP - 0 



Integer designating the vertical axes of 
the viewing plane where 1 = X<j, 2=Y^,and 3=Zo . 
(see Figure 5.) 

Angular rotation of model about its X-axis in 
degrees (must be performed third, see Figure 

5.) 

Angular rotation of model about its Y-axis in 
degrees (must be performed second, see Figure 
5.) 

Angular rotation of model about its Z-axis in 
degrees (must be performed first, see Figure 
5.) 

1 Frame change before plotting 
(a frame change resets the Y-origin past 
previous plot by XSPACE given in Namelist 
OPTION and resets the X-origin at 0.0) 

0 no frame change before plotting 

1 automatic scaling of plot and computation 
of proper origin location 

2 user-specified origin and scaling 
Maximum dimension desired on completed plot, 
in inches (used for scaling if ISCALE=1) 
X-location of plot origin 

(used if ISCALE=2) 

Y-location of plot origin 
(used if ISCALE=2) 

Model size reduction factor 

(i.e., PSCALE is equal to actual model size 

divided by desired plot size, used 

if ISCALE=2) 

0 no numbering on plots 

1 numbering of grid points 

2 numbering of elements 

Height of integers specified by NOTAT, in 
inches (must be > 0.07) 

0 plot of undeformed structure 
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IDMAG - 2 



DMAGS - 1.0 



KSYMXY - 0 
KSYMXZ - 0 
KSYMYZ - 0 



1 plot of deformed structure 

2 exploded plot 

3 displacements represented by vectors 

1 direct magnification of displacement data 
by DMAGS 

2 scaling of displacement data to a 
maximum value ofDMAGS 
Magnification of displacements 

(if KDISP= 1 or 3) 

Reduction factor of elements 
(if KDISP=2) 

1 symmetry about X-Y plane 

1 symmetry about X-Z plane 

1 symmetry about Y-Z plane 



Symmetries are performed consecutively {±.e., a plate 
quadrant with KSYMXZ and KSYMYZ equal to one would yield a 
complete plate ) . 



XXMAX,XXMIN - 1.0E+20 Locate cutting planes parallel to 
principal planes 

YYMAX,YYMIN - 1.0E-20 (X-Y, X-Z, Y-Z) to limit plot 

ZZMAX,ZZMIN 



NDMAX - 9999999 
NDMIN - 0 
NELMAX - 9999999 
NELMIN - 0 
KODE - 0 



Maximum grid point identification 
number to be included in the plot 
Minimum grid point identification 
number to be included in the plot 
Maximum element identification 
number to be included in the plot 
Minimum element identification 
number to be included in the plot 
Specifies control option after a 
plot is complete 

0 last plot, exit from program 

1 read another Namelist PICT 
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2 read a new set of displacement 
data (see NOTE 1 below) 

3 read a complete new set of input 
data starting with a title card. 

(NOTE 1: For SAP IV displacement data, K0DE=2 signifies 
that the next load case displacements will be assigned to 
the model nodal points.) 

The previous sections describe a complete basic set 
of input data, if KODE =0 in Namelist PICT. For KODE =1, 2, 
or 3, additional sections of the deck must be repeated. The 
deck must end with a Namelist PICT having a value of KODE = 
0 in it. An example input data deck and output plors for 
the simple truss problem of Appendix A is found in Appendix 
B. 



C. METHOD FOB ALTERING SUBROUTINE PSAP 



In the event a user has geometry and displacement data 
decks from a program other than SAP IV, it is possible to 
plot those decks with PSAP. The subroutine will handle 
rod-like elements, triangular elements, or quadrilateral 
elements when they are input in acceptable format. By 
studying subroutines GE0M9 and DATA9, and SUBROUTINE PSAP 
(Appendix D) , the necessary sequence of input can be 
determined. PSAP is presently constructed so that a user 
may supply his own routines through the use of subroutines 
GE0M1 or GEOM2 for geometry data and DATAI or DATA5 for 
displacement data. In order to add a subroutine to PSAP, it 
need only be placed after the main calling program in the 
sequence of control cards as discussed in paragraph II. B. 2. 
The essential features for the input of the geometry data 
are the nodal points, with their X,Y,Z coordinates, and the 
connectivity sequence for the finite element model. The 
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necessary part of displacement data input is the node point 
number, with the 0,V,M displacements. Adding a user-prepared 
subroutine through GE0M1 or 2 and DATAl or 5 should prove to 
be a relatively straight-forward taslc for the user who 
desires to do so. The listing of GE0H9 in Appendix D can be 
used as a guide. 



D. SIGNIPICANT ASPECTS OF SUBEOOTINE PSAP 



1 • Exploded Plots 

Often the absence or presence of elements in a 
finite element model cannot be determined from a 
conventional oblique orthographic projection. For example, a 
line element that is coincident with an edge of a triangular 
or quadrilateral element could not be detected. To show 
clearly each element, PSAP contains an algorithm for 
generating exploded oblique orthographic projections. This 
can be a valuable tool in checking the toplogy of an 
analytical model, (i.e., KDISP = 2) 

2- Portions of Mode ls 

The ability to isolate a portion of a model for detailed 
examination is a very useful and desirable asset of a 
preprocessor. SOBEOUTIME PSAP has the capability of 
specifying cutting planes (i.e.,XXMIN ,XXMAX , YYMIN , YYMAX , 
etc.) or maximum and minumum element numbers (i.e., NELMIN, 
NELMAX) . Examples of this are shown in Appendix D. 
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3. 



Specifi cation of View 



The specif ication of view of a model is done through 
the use of the parameters KVERT,KHORZ, PSI, THETA, and PHI 
that are found in Namelist PICT (Section III.B.9) The 
specific details of how this is accomplished can be found in 
Reference 2. There are a great number of possible 
combinations of the above parameters, and Appendix C 
illustrates various combinations that were used on a section 
of a finite element wing model. The Calcomp plots, along 
with the parameters as specified, are found in each figure. 



^ • P ossible Sources of E rr ors 



(а) The most probable source of error could be 
incorrect data deck preparation or deck sequencing. It is 
important to use Figure 2 as a guide while preparing the 
input for PSAP. 

(б) Other errors may occur if the arrays ZZZ and 
DISPD are not dimensioned correctly as discussed in 
paragraph III.B.3. During the execution of the program 
several manipulations are performed with the two arrays, and 
it is possible for addresses of the data to be lost within 
the IBM 360/67. The 360 System error messages may not 
directly indicate that the dimensions of the arrays have 
been exceeded. 

(c) A most important point to remember in generating 
a sequence of plots is that once a parameter has been 
assigned a value in a namelist, it retains that value until 
it is reassigned. For example, if PLOTSZ is assigned a value 
of 8.0 for the first in a series of plots and it is not 
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redefined in any subsequent Namelists PICT , the value of 
PLOTSZ will retain the value (8.0), as originally specified. 

(d) It is possible to make any number of errors, 
however, and all of them cannot be anticipated. The error 
messages from the IBM 360/67 System will, in most cases, be 
straight-forward and facilitate easy trouble-shooting. 
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IV. COHCLOSIOHS AND REC OMM END ATIONS 



General purpose finite element structural analysis 
programs like SAP IV are significant analytical tools that 
can be used most effectively when coupled with a flexible 
pre- and postprocessor. The availability of a powerful 
structural analysis tool with partial pre- and 
postprocessing capability has now been provided for Naval 
Postgraduate School Students. However, at present this 
capability is limited to pre- and postprocessing the data of 
only five of the SAP IV elements. It does not appear 
feasible, at this time, to add any more capability as far as 
the number of SAP IV elements goes. However, other 
additional improvements to postprocessing at NPS should be 
made. For example. Appendix C of Reference 2 lists a 
program developed to produce contour plots of stress data 
from finite element models. The program, in all probability, 
could be adapted to SAP IV and the NPS IBM 360/67. This 
represents a possible avenue to further the present 
structural analysis capability for NPS students. SAP IV 
itself can be expanded. There is a great deal of work being 
done in the area of composite materials, and SAP IV 
possesses only very limited orthotropic material capability. 
However, there are subroutines within SAP IV that could be 
modified to reflect more current methods in handling 
composite materials. 

In conclusion, the finite element method in structural 
engineering is the most powerful analysis tool available 
today and needs to be exploited by NPS students in both 
classroom and research work. Programs like SAP IV are 
available at many Navy laboratories, development centers. 
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and aerospace companies throughout the country. The current 
trends in structural analysis are toward computer-aided 
techniques, and exposure to a general purpose package for 
students enrolled in AE 4101 and AE 4102 is a significant 
asset . 
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APPENDIX A 



EXAMPLE PROBLEM USING SAP IV 



The schematic given below illustrates the example 
problem for which the data on the following page were 
prepared . 




E = 10.1E+06 psi 
ALPHA = 12. 6E-06 psi 
A = 0.5 sq.in. 

A = 1.0 sq.in. 

\ = 0.75 sq.in. 
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STANCARC NPS EOF CARD 



APPENDIX B 



EXAMPLE PROBLEM USING SUBROUTINE PSAP 



The pin-jointed truss problem in Appendix A is also used 
as an example in illustrating SUBROUTINE PSAP use. The 
entire computer card deck used to generate Figures 6, 1 , and 
8 is listed on the following two pages of this appendix. 
The printed computer output for this problem is shown in 
Figure 9. Note that the displacements of node number 4 have 
been magnified by 1000. 
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Figure 6 - UNDEFORMED TRUSS MODEL 
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Figure 7 - DEFORMED TRUSS MODEL<NODE4 DISPLACED^ 
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Figure 8 - DISPLACEMENTS < N 0DE4 > FOR TRUSS MODEL PLOTTED 
AS A VECTOR 
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Figure 9 - SAMPLE COMPOTER PRINTOUT FROM PSAP OF TRUSS 
MODEL 
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APPENDIX C 



SAMPLE SUBROUTINE PSAP OUTPUTS FOR A FINITE ELEMENT MODEL OF 
A WING 



The figures contained in this appendix are of a finite 
element wing model developed in . Flight Vehicle Structural 
Analysis II (AE 4102, QTR.IV, 1976). The original wing was 
designed in AE 4274, which is the latter half of the 

subsonic structural design sequence in the Department of 
Aeronautics. The finite element model of the wing was 
developed in AE 4102 as a class project. (The structural 
analysis and design sequences in the the Department of 
Aeronautics are constructed so that they interelate the two 
different areas.) Each student was assigned responsibility 
for a portion of the wing and the student whose design was 
used was the analysis leader. 

Figure 10 shows the axis orientation used in the 

original definition of the wing as designed in AE 4274. 

Figure 11 and 12 show the entire model with over four 

hundred membrane elements and with different views of the 
model specified. In Figures 13 and 14 a single portion of 
the wing is shown with different rotation angles. The 
significant Namelist PICT values used to orient the model 
are included in the figures. 
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Figure 10 - WING HODEL AXIS ORIENTATION 
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KH0RZ=1 , KVERT=3 , PHI=70 
THETA=-25,PS I=-10 




Figure 11 - ORTHOGRAPHIC PROJECTION OF WING MODEL -1 
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KH0RZ=2 , KVERT=3 , PHI=60 
THETA=-25,PSI=-10 




Figure 12 - 



ORTHOGRAPHIC PROJECTION OF WING MODEL -2 




Figure 13 - wing MODEL SECTION - ORIENTATION -1 
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APPENDIX D 



SUBROUTINE PSAP LISTING 
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